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We investigate the problem of speciation and coexistence in simple ecosys- 
tems when the competition among individuals is included in the Eigen model 
for quasi-species. By suggesting an analogy between the competition among 
strains and the diffusion of a chemical inhibitor in a reaction-diffusion system, 
the speciation phenomenon is considered the analogous of chemical pattern 
formation in genetic space. In the limit of vanishing mutation rate we ob- 
tain analytically the conditions for speciation. Using different forms of the 
competition interaction we show that the speciation is absent for the genetic 
equivalent of a normal diffusing inhibitor, and is present for shorter-range 
interactions. The comparison with numerical simulations is very good. 
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In this work we address the problem of speciation (species formation) in simple ecosys- 
tems, mirroring some aspects of bacterial and viral evolution. Our model can be considered 
as an extension of the Eigen model With respect to the latter, we introduce the 

competition among individuals. 

Eigen's phenomenological theory of self-reproducing macromolecules (or haploid organ- 
isms) illustrates the concept of stable quasi-species, i.e. a peaked distribution of genomes 
around a master sequence, its width being determined by mutations. In its simpler for- 
mulation, the various genomes have different reproductive rates, the logarithm of which 
constitutes the fitness landscape The master sequence is located in correspondence 

of the maximum of the fitness. In general a one to one correspondence between a given 
phenotype and a genotype is assumed (no polymorphism nor age structure). The genomes 
are coupled by mutations and by a global constraint on the total number of individuals 
(constant organization). One usually considers only point mutations (the most common 
ones), which correspond to a diffusion process in genetic space. In this way one can define 
the concept of distance in genetic space as the number of mutations needed to connect two 
genomes. The Eigen model has also been studied in the contest of statistical mechanics 



For a vanishing mutation rate (which is the case for actual organisms), only one quasi- 
species can survive in the long time limit, except for the marginal case of degenerate maxima 
of the fitness or for an extremely rough landscape (similar to the spin glass energy landscape) 
for which the system never attains equilibrium. Epstein studied the problem without 



considering mutations; he showed that the coexistence is possible if the species are self- 
limiting (i.e. there exists a form of self-competition, modeled for instance by a logistic term) 
and coexisting species does not compete directly. On the contrary, when two species are in 
competition (because they share some resource — an enzyme in Epstein's case), only the 
fittest one survives. However, he did not introduce the genetic distance among species nor 
presented any evolutionary mechanism for the speciation phenomenon. 

We think that the direct competition for local resources among strains, coupled with 



a weak mutation rate, is the simplest mechanism for modehng both speciation and stable 
coexistence in simple smooth landscapes. The mutations are needed to populate newly 
formed niches, while the competition actively separates the strains into quasi-species. One 
can consider the following analogy with a Turing mechanism for chemical pattern formation. 
The main ingredients are an autocatalytic reaction process (reproduction) with slow diffusion 
(mutations) coupled with the emission of a short-lived, fast-diffusing inhibitor (competition). 
In this way a local high concentration of autocatalytic reactants inhibits the growth in its 
neighborhood, acting local negative interaction. 

In genetic space, the local coupling is given by the competition among genetically kin 
individuals. For instance, assuming a certain distribution of some resource (such as some 
essential metabolic component for a bacterial population), then the more genetically similar 
two individuals are, the wider the fraction of shared resources is. The effects of competition 
on strain x by strain y are modeled by a term proportional to the relative abundance of the 
latter, p{y), modulated by a function that decreases with the genetic distance between x 
and y. Another example of this kind of competition can be found in the immune response in 
mammals. Since the immune response has a certain degree of specificity, a viral strain x can 
suffer from the response triggered by strain y if they are sufficiently near in an appropriate 
genetic subspace. Again, one can think that this effective competition can be modeled by 
a term, proportional to the relative abundance of the strain that originated the response, 
which decreases with the genetic distance. 

Although Eigen's model is defined in an hypercubic genetic space, and the error threshold 
transition rigorously exists only in an infinite-dimensional space , the proposed speciation 
phenomenon is independent on the dimension of the genetic space. We shall work therefore 
in a linear genetic space. An instance of a similar (sub-) space in real organisms is given by 
a repeated gene (say a tRNA gene): a fraction of its copies can mutate, linearly varying the 
fitness of the individual with the "chemical composition" of the gene This degenerate 
case has been widely studied (see for instance Refs. |I2[); one should introduce the multi- 
plicity of a degenerate state, which can be approximated to a Gaussian, but if one works in 



the neighborhood of its maximum (the most common chemical composition) the multiphcity 
factors are nearly constants. Another example is given by the level of catalytic activity of 
a protein. A linear space has also been used for modeling the evolution of RNA viruses on 



HeLa cultures [13 



Let us start with a one dimensional "chemical" model of cells that reproduce asexually 
and slowly diffuse (in real space), p = p{x,t) being their relative abundance at position x 
and at time t. These cells const it utively emit a short-lived, fast- diffusing mitosys inhibitor 
q = q{x, t). This inhibitor may be simply identified with some waste or with the consumption 
of a local resource (say oxygen). The diffusion of the inhibitor is modeled as 

Qq d'^q 

^ = ^ol> + I^^-M, (1) 

where ko, ki and D are the production, annihilation and diffusion rates of q. 
The evolution of the distribution p is given by 

| = (^(^'^)-^W)p + /^0' (2) 

~A{t) = J A{y,t)p{y,t)dy. (3) 

The growth rate A can be expressed in terms of the fitness H as 

A{x,t) = exp{H{x,t)) . (4) 

Due to the form of equation (^, the distribution p is always normalized to one. The diffusion 
rate of q, D, is assumed to be much larger than fi. The growth rate A, can be decomposed 
in two factors, A{x,t) = AQ{x)Ai{q{x,t)), where Aq gives the reproductive rate in absence 
of q, so Ai(0) = 1. In presence of a large concentration of the inhibitor q the reproduction 
stops, so Ai(oo) = 0. A possible choice is 

A{x, t) = exp{Ho{x) — q{x, t)). 

For instance, Hq{x) could model the sources of food or, for algae culture, the distribution 
of light. 



Since we assumed a strong separation in time scales, we look for a stationary distribution 
q{x,t) of the inhibitor (eq. (|I])) by keeping p fixed. This is given by a convolution of the 
distribution p: 



q{x,t) = J Jexp (- ^^ ^^^ jp{y,t)dy, 



where J and R depend on the parameters k^, ki, D. In the following we shall use J and R 
as control parameters, disregarding their origin. 

We can generalize this scenario to non-linear diffusion processes of the inhibitor by using 
the react ion- diffusion equation eq. (^, with the fitness H and the kernel K are given by 

H{x, t) = Ho{x) -JJk [^^) P{y, t)dy (5) 

K{r) = exp , (6) 

i.e. a symmetric decreasing function of r with K{0) = 1. The parameters J and a control 
the intensity of the competition and the steepness of the interaction, respectively. 

Let us consider the correspondence with the genetic space: the quantity x now identifies 
a genome, the diffusion rate n is given by mutations, and the inhibitor q (which is no more 
a real substance) represents the competition among genetically related strains. The effects 
of competition are much faster than the genetic drift (mutations), so that the previous 
hypotheses are valid. While the genetic interaction kernel K{r) is not given by a diffusion 
process, its general form should be similar to that of eq. @: a decreasing function of the 
genetic distance between two strains. We shall refer to the p-independent contribution to 
the fitness, Hq{ the static fitness landscape. 

Our model is thus defined by eqs. (0-|^). We are interested in its asymptotic behavior 
in the limit fi ^ 0. Actually, the mutation mechanism is needed only to define the ge- 
netic distance and to allow population of an eventual niche. The results should not change 
qualitatively if one includes more realistic mutation mechanisms. 

Let us first examine the behavior of eq. (^ in absence of competition (J = 0) for a 
smooth static landscape and a vanishing mutation rate. This corresponds to the Eigen 



model in one dimension: since it does not exhibit any phase transition, the asymptotic 
distribution is unique. The asymptotic distribution is given by one delta function peaked 
around the global maximum of the static landscape, or more delta functions (coexistence) 
if the global maxima are degenerate. The effect of a small mutation rate is simply that of 
broadening the distribution from a delta peak to a bell-shaped curve . 



While the degeneracy of maxima of the static fitness landscape is a very particular 
condition, we shall show in the following that in presence of competition this is a generic 
case. For illustration, we report in Figure 1 the numerical computation of the asymptotic 
behavior of the model for a possible evolutive scenario that leads to the coexistence of 
three species. We have chosen a smooth static fitness Hq (see eq. (^) and a Gaussian 
{a = 2) competition kernel. The effective fitness H is almost degenerate (here fi > and 
the competition effect extends on the neighborhood of the maxima), and this leads to the 
coexistence. One could show that the curvature of the maxima affects the width and the 
height of the quasi-species distribution in presence of mutations Jl^ . 



We shall now derive the conditions for the coexistence of multiple species. Let us assume 
that the asymptotic distribution is formed by L delta peaks pk, k = 0, ... ,L — 1, for a 
vanishing mutation rate (or L non-overlapping bell shaped curves for a small mutation rate) 
centered at i/k- The weight of each quasi species is jk, i-e. 



„ L-l 

/ pk{x)dx = -fk, 



k=0 

The quasi-species are ordered such as 70 > 71, . . . , > 7l-i- 
The evolution equations for the pk are (/x 0) 

dpk 



dt 



{A{yk) - A)pk, 



where A[x) = exp {H{x)) and 



- y, 



H{x) = Hoix) 'JUK [-^) Ir 

The stability condition of the asymptotic distribution is {Aiijk) — A)pk = 0, i.e. either 
A[yk) = A = const (degeneracy of maxima) or = (all other points). In other terms one 



can say that in a stable environment the fitness of all individuals is the same, independently 
on the species. 

The position yk and the weight 7*; of the quasi-species are given hy A{yk) = A = const 
and dA{x) / dx\ — 0, or, in terms of the fitness H, by 



Let us compute the phase boundary for coexistence of three species for two kinds of 
kernels: the exponential (diffusion) one {a = 1) and a Gaussian one (a = 2). 

We assume that the static fitness Hq{x) is a symmetric linear decreasing function except 
in the vicinity oi x — 0, where it has a quadratic maximum: 



so that close to x — one has Ho{x) ~ —hx'^/r^ and for x ^ 00, Hq{x) ~ 6(1 — \x\/r). We 
have checked numerically that the results are qualitatively independent on the exact form 
of the static fitness, providing that it is a smooth decreasing function. 

Due to the symmetries of the problem, we have one quasi-species at x = and two 
symmetric quasi-species at a; = ±|/. Neglecting the mutual influence of the two marginal 
quasi-species, and considering that H'q{Q) = K'{Q) = 0, K'{y/R) = -K'{-y/r), K{0) = J 
and that the three-species threshold is given by 70 = 1 and 71 = 0, we have 



where y — y/R, f — r/R and b — b/ J. In the following we drop the tildes for convenience. 






(7) 




^ + K'{y) = 0. 



Thus 
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r — z — G exp 




Gz'^-'exp (-—) = 1, 



V a J 



where G = r/b. 



For a = 1 we have the coexistence condition 



ln(G) = r-l + G. 



The only parameters that satisfy these equations are G = 1 and r = 0, i.e. a flat landscape 
(6 = 0) with infinite range interaction {R = oo). Since the coexistence region reduces to a 
single point, it is suggested that a = 1 is a marginal case. 
For a = 2 the coexistence condition is given by 



One can solve numerically this system and obtain the boundary Gc{r) for the coexistence. 
In the limit r — > (static fitness almost fiat) one has 



with G,(0) = 2.216. . .. 

We have performed several numerical simulations for different values of the parameters, 
whose results are presented in Figure 2. The boundary of the multi-species phase is well 



rate /i, at least for /i < 0.1, which can be considered a very high mutation rate for real 
organisms. The most important effect of /i is the broadening of quasi-species curves, which 
can eventually merge. 

In conclusion, we have introduced a model for the genetic evolution of haploid organisms 
under the pressure of a static fitness landscape and competition. This model exhibits the 



-{G + r)z + 1 



0, 




G,(r) ^ G,(0) - r 



(8) 



approximated by eq. (H); in particular, this boundary does not depends on the mutation 
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phenomenon of species formation in a way reminiscent of a chemical pattern formation via a 
Turing-hke mechanism. We have analyzed analytically this system in the limit of vanishing 
mutation rate and linear genetic space, showing that an increasing level of a short-range 
competition induces a transition from a single species distribution to a stable environment 
in which multiple genetically distinct species are present. The comparison of the analytical 
approximation with the numerical integration of the original differential equations is very 
good. We think that the mechanism that we proposed is the simplest one for modeling 
speciation and species coexistence in a smooth (or fiat) fitness landscape. 

We wish to thank G. Guasti, G. Cocho, L. Peliti, G. Martinez-Mekler and P.Lio for fruit- 
ful discussions. M.B. thanks the Dipartimento di Matematica Applicata "G. Sansone" for 
friendly hospitality. Part of this work was done during the workshop Chaos and Complexity 
at ISI- Villa Guahno (Torino, Italy) under CE contract ERBCHBGCT930295. 
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FIGURE CAPTIONS 



Fig. 1. Static fitness Hq, effective fitness H, and asymptotic distribution p numerically 
computed for the following values of parameters: a = 2, /i = 0.01, Hq = 1.0, b = 0.04, 
J = 0.6, R = 10, r = 3 and N = 100. 

Fig. 2. Three-species coexistence boundary Gc for a = 2. The continuous line repre- 
sents the analytical approximation, eq. (^, the circles are obtained from numerical 
simulations. The error bars represent the maximum error. 
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